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SUMMARY 

The purpose of this paper is to develop the Dugdale model for the compact 
specimen. In particular, plastic-zone size and crack-tip-opening displacement 
(CTOD) equations were developed. Boundary-collocation analyses were used to 
analyze the compact specimen subjected to various loading conditions (pin 
loads, concentrated forces, and uniform pressure acting on the crack surface). 
Stress-intensity factor and crack -surface displacement equations for some of 
these loadings were developed and used to obtain the Dugdale model. The 
results from the equations for plastic-zone size and CTOD agreed well with 
numerical values calculated by Terada for crack-1 ength-to-width ratios greater 
than 0.4. 


Associate Professor, University of Missouri, Rolla, MO 65401. 
Senior Scientist, NASA Langley Research Center, Hampton, VA 23665. 





2 


NOMENCLATURE 

a crack length measured from centerline of holes 

Bj ratio of bj to specimen width (bj/w) 

b distance from crack tip to concentrated force 

bj dimensions for partially-loaded crack (j = 1,2) 

c distance measured from hole centerline to concentrated force 

d crack length plus tensile plastic zone 

E Young's modulus of elasticity 

F boundary-correction factor on stress intensity for pin-loaded hole 

6 boundary-correction factor on stress intensity for concentrated force 

H boundary-correction factor on stress intensity for uniform pressure 

K stress-intensity factor 

P pin load per unit thickness on holes 

Q concentrated force per unit thickness on crack surface 

r radius of circular holes 

V crack-surface one-half displacement 

w specimen width 

w‘ specimen width measured from location of concentrated force 

x,y Cartesian coordinates 

z complex variable, z = x + iy 

a crack-1 ength-to-width ratio (a/w) 

B plastic-zone-size-to-width ratio (p/w) 

Y equal p/(a + p) 

A equal b/w' or b/(w - a + b) 

n material constant, n = 0 for plane stress and n = v for plane 

strain 


X 


equal (a + p)/w 
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v Poisson's ratio 

£ distance measured along crack surface from crack tip 

p length of tensile plastic zone 

a uniform pressure acting on segment of crack surface 

o 0 flow stress of material 

$ ,¥ complex stress functions 


INTRODUCTION 

Many aspects of material behavior contribute to the growth of cracks in 
plates under monotonic and cyclic loading. A general view concerning the 
behavior of material at the crack tip is that the effects of plastic flow must 
be accounted for to correlate and to predict the growth of cracks. The Dugdale 
model [1] is a very simple approach that simulates the effects of plastic flow 
on plastic-zone size and on crack-tip-opening displacements (CTOD). 

The Dugdale model concept has been used with center-crack tension speci- 
mens in several fracture analyses [2-4] and in the development of fatigue 
crack-closure models [5,6] for cracked metallic materials. Currently, the com- 
pact specimen is the most widely used specimen for fatigue crack growth rate 
and fracture toughness testing. The purpose of this paper is to develop the 
Dugdale model for the compact specimen (Fig. 1(a)). In particular, plastic- 
zone size (p) and CTOD equations have been developed. The Dugdale model is the 
superposition of two elastic problems. These two elastic problems, a compact 
specimen subjected to pin loads (P) or to uniform pressure (o 0 ) applied to a 
segment of the crack surfaces, are shown in Fig. 1(b). To develop the Dugdale 
model for the compact specimen, additional equations for stress-intensity 
factor (K) and crack-surface displacement (V) must be obtained for this speci- 
men under various loading conditions. These loadings are shown in Fig. 2. 



The pin-loaded case. Fig, 2(a), was previously analyzed by Newman [7] to 
obtain stress-intensity factors and crack-surface displacements. Srawley [8] 
fitted the results from Ref 7 with a power series to obtain an equation for K. 
The boundary-collocation method [7] was used herein to analyze a pair of con- 
centrated forces on the crack surfaces (Fig. 2(b)) and a uniform pressure 
acting on a segment of the crack surfaces (Fig. 2(c)). Both stress-intensity 
factors and crack-surface displacements under plane-stress conditions (v = 0.3) 
were calculated. Stress -intensity factor equations were developed for these 
loadings. Crack-surface displacement equations were also developed for the 
pin-loaded case (Fig. 2(a)) and for the uniform-pressure case (Fig. 2(c)) with 
b 2 = 0. 

The stress-intensity factor equations for pin loading and uniform pressure 
(b 2 = 0) were used to develop the Dugdale model (plastic-zone size) for the 
compact specimen. Again, an equation was fitted to these results which gives 
the plastic-zone size as a function of crack length, specimen width, applied 
load, and flow stress. Likewise, an equation was developed for the crack-tip- 
opening displacement (CTOD) for the Dugdale model. The results from these 
equations were compared with recent numerical results derived by Terada [9]. 

ANALYSIS OF THE COMPACT SPECIMEN 

The boundary-collocation method [7,10] was used to analyze the compact 
specimen subjected to a pair of concentrated forces, as shown in Fig. 2(b). 

Some details of the analysis are given in the Appendix. The stress-intensity 
factors for a wide range of configuration parameters (a/w and b/w') were cal- 
culated. The following empirical equation was then fitted to the numerical 
results 
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where a = a/w and A = b/w'. The correction factor, G, will be presented 
later. 

The crack-surface displacements for the compact specimen subjected to the 
three types of loading shown in Fig. 2 were calculated by a procedure suggested 
by Tada, Paris, and Irwin [11]. This procedure is summarized here. 

Consider a cracked body loaded by a force system, p, as shown in Fig. 3. 
Let a pair of virtual forces, q, be applied on the crack surface at the loca- 
tion where the crack-surface displacement is desired. If K p and K q are 
stress-intensity factors due to p and q, respectively, then the crack- 
surface displacement, V, at location x = c - d is 


1 n 2 d 9K n 


( 2 ) 


where n = 0 for plane stress and n = v for plane strain. For the compact 
specimen, the partial derivative, 3K q /3q, was obtained from Eq (1) as 

_ 2G(ct,A) _ 2G(ot,A) 

3Q /5b " /2i(a - c) C 

where b is replaced by (a - c) in the function G. During the integration, 
the crack length c in Eq (2) is held constant (see Fig. 3) while the crack 
length (a) changes from c to d. The crack length d is the final crack 
length for which the displacement is desired. Because K p and K q are 
generally complicated functions of crack length, Eq (2) was numerically 
integrated. 


4 
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STRESS-INTENSITY FACTORS 
Pin-Loaded Holes 

As previously mentioned, the standard compact specimen subjected to pin 
loading, as shown in Fig. 2(a), has been previously analyzed [7], Stress- 
intensity factors were calculated for crack-1 ength-to-width ratios ranging from 
0.2 to 0.8. Srawley [8] then fitted a polynomial expression to these results 
and it was 

K P=^ F (t) <«> 

where 

F(ct) = (2 + a) (0.886 + 4.64a - 13.32a 2 + 14.72a 3 

- 5.6a 4 )/ ( 1 - a) 3/2 (5) 

and a = a/w. Equation (5) is within ±0.3 percent of the collocation results 
for 0.2 a < 1 . Equation (4) also approaches the exact asymptotic solution 
as a approaches unity [12]. 


Concentrated Forces 

The standard compact specimen subjected to a pair of concentrated forces 
(see Fig. 2(b)) was analyzed by the boundary-collocation method. The stress 
functions used to exactly satisfy the crack-surface boundary conditions and to 
approximately satisfy the boundary conditions on the holes and the external 
boundary are discussed in the Appendix. Stress-intensity factors were calcu- 
lated for a wide range of configuration parameters. An empirical equation was 
then fitted to the numerical results and is 



( 6 ) 
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where 

G(a,A) = (l + A x A + A 2 A 2 )/(l - A ) 3 ^ 2 
A x = 3.57 + 12.5(1 - a) 8 
A o = 5.1 - 15.32a + 16.58a 2 - 5.97a 3 


(7) 

(8) 

(9) 


a = a/w, A = b/w' = b/(w - a + b), and 0 £ a < a. Equation (6) is within 
±0.5 percent of the collocation results for 0.3 £ a £ 0.8 and within ±1 per 
cent for 0.2 < a < 0.3. 


Uniform Pressure Acting on a Segment of Crack Surface 
With the availability of the K-solution for concentrated forces, the 
stress-intensity factors for the standard compact specimen subjected to a uni- 
form pressure, o, acting on a segment (b^ to bg) of the crack surface, 
shown in Fig. 2(c), can be determined using the Green's function concept. The 
stress-intensity factor is then given by 




( 10 ) 


Integration of Eq (10) for uniform pressure gives 



(ID 


H = 


(1 


— 3/2 [2B(l + Aj + A 2 ) \|b 2 + (1 - a)B 


+ (1 - a)(5 + Aj -3A 2 )\|b 2 + (1 - a)B 

+ (1 - a) 2 (3 - Aj + 3A 2 ) An + /B + 1 - a)l| B Bl 

J IB = Bo 


( 12 ) 
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where a = a/w and Bj = bj/w. The functions Aj and A2 are given 
by Eqs (8) and (9), respectively. Because the concentrated force solution was 
within ±0.5 percent of the collocation results, the uniform pressure case would 
have the same accuracy for 0 £ Bj £ a and 0.3 £ a _< 0.8. Equation (11) is 
useful in developing the Dugdale model and plastic-zone size equations for the 
compact specimen. This development will be discussed later. 

CRACK-SURFACE DISPLACEMENTS 


Pin-Loaded Holes 

The crack-surface displacements for the compact specimen under pin loading 
(Fig. 2(a)) for plane-stress conditions were previously calculated using a 
collocation analysis [7,13], For comparison, the displacements were also cal- 
culated from Eq (2). The crack-surface displacement, V p , under plane-stress 
conditions at a distance c from the hole centerline is 


V -|E. 7 F(«)6(«.A) .. 

^ c/w /2 tt(o - c/w) 


(13) 


Equation (13) was numerically integrated. Because the polynomial expressions 
for F and G (Eqs (5) and (7), respectively) are within ±0.5 percent of the 
collocation results for 0.3 _< a 0.8, then the computation of crack-surface 
displacement would, in general, be within ±1 percent of the collocation 
results. The normalized displacements, EV p /P, as functions of 5/a, are 
shown in Fig. 4 for various values of a/w. The synt>ols show calculations from 
Eq (13) and the curves were obtained from a boundary-collocation analysis [13]. 
The results from numerical integration and collocation agreed within 1 percent. 
Equation (13) cannot be used for locations 5/a > 1 - 0.2 (w/a) because F and 
G were not developed for a/w ratios less than 0.2. 




i 


% 
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For ease of computation, a polynomial expression for the normalized dis- 
placements shown in Fig. 4 was developed. This expression is 


f!e. = 4 |_l 

p “l o_. 


1 +C, i + C*± 


2irw L A ' nva J 2\a/ J 


C x = -1.25 + 9.76a - 20.15a 2 + 16.62a 3 


C 2 = 0.64 - 4.34a + 10.24a 2 - 8.46a 3 


(14) 


where a = a/w and F is given by Eq (5). Equation (14) is within ±1 percent 
of the numerical results from Eq (13) for 5/a < 1 - 0.2/a and 0.3 < a < 0.8. 


Concentrated Forces 

The crack-surface displacements for the compact specimen subjected to a 
pair of concentrated forces (Q) (Fig. 2(b)) under plane-stress conditions were 
calculated using either the collocation analysis (see Appendix) or Eq (2). 
Using Eq (2), the crack-surface displacement, Vq, at a distance c from the 
hole centerline is 


v = 2Q « 6(°.A 0 )6(»,t) 

^ c/w J(a - c/w)b Q /w 


(15) 


where A 0 = b 0 /(w - a + b 0 ) in Eq (7). The distance b 0 is the loca- 
tion of the stationary force Q. The distance b is the location of the 
virtual force q. Again, b must be replaced by (a - c) and b 0 must be 
replaced by (a - c 0 ). The distance c 0 is the distance from the hole 
centerline to the location of the force Q. Again, Eq (15) was numerically 
integrated. 

The normalized displacements, EVq/Q, as functions of 5/a, are shown in 
Fig. 5 for various values of a/w. The symbols show calculations from Eq (15), 
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and the curves were obtained from the boundary-collocation analysis. The 
results from numerical integration and collocation again agreed well. Again, 
Eq (15) cannot be used for locations 5/a > 1 - 0.2(w/a) because G was not 
developed for a/w ratios less than 0.2. 


Uniform Pressure Acting on a Segment of Crack Surface 
The crack-surface displacements for the compact specimen subjected to a 
uniform pressure, as shown in Fig. 2(c), were also calculated using Eq (2). 
The displacement, V a , at a distance c from the hole centerline is 


V 


a 


ow a r H(“.B 1 ,B 2 )G(a,A) 

orr J , ■ - ■ a<* 

c/w /a - c/w 


(16) 


where G and H are given by Eqs (7) and (12), respectively. The function 
Bj is bj/w, and bj must be replaced by (a - Cj). The distance 
cj is measured from the hole centerline to appropriate ends of the uniform 
pressure loading. 

Some typical normalized crack-surface displacements, EV 0 /(c w), with 
a/w = 0.6 for various lengths of uniform pressure are shown in Fig. 6. The 
displacements calculated from Eq (16), shown as symbols, were within 1 percent 
of displacements calculated from collocation (curves) because the polynomial 
expressions for G and H (Eqs (7) and (12), respectively) are within 
±0.5 percent of collocation results. These displacements are useful in devel- 
oping the CTOD equations for the Dugdale model. 




DUGDALE MODEL FOR COMPACT SPECIMEN 

The Dugdale model for the compact specimen (Fig. 1) requires that the 
"finiteness" condition of Dugdale be satisfied. This condition states that the 
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K at the tip of the plastic zone (at d = a + p) is zero. From this condi- 
tion, the plastic-zone size (p) was calculated for various a/w and P/(wa 0 ) 
ratios. An equation was then fitted to these results and was 

p _ uw/PF \ 2 

a 8a(wo o j o (17) 

where 

F = (2 + a) (0.886 + 4.64a - 13.32a 2 + 14.72a 3 

- 5.6a 4 )/(l - «)3/2 

fj = -2.687(1 - a) + 0.167/a 2 
f 2 = -2.48 - 0.039/(1 - a) 6 

and a = a/w. Equation (17) is within ±1 percent of the current collocation 
results for 0.3 < (a + p)/w < 0.8 and p/(w - a) < 0.5. A comparison of 
plastic-zone sizes calculated from Eq (17) and some recent collocation results 
from Terada [9] are shown in Fig. 7. The normalized plastic-zone size, 
p/(w - a), is plotted against P/(wa 0 ) for various values of a/w. The 
curves show the present results (Eq (17)) and they agreed well with the previ- 
ous results from Terada (symbols). 

In the compact specimen, the material at point A in Fig. 1(a) is in com- 
pression. At a certain load, this material will yield in compression. From a 
finite-element-strip-yield analysis for an elastic-perfectly plastic material. 
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the load that causes incipient yielding at point A was calculated. For each 
value of a/w, these loads are shown as the upper ends of the solid curves. 

The upper ends of the dashed curves denote loads and plastic-zone sizes where 
backface yielding at point A has about a 2 percent effect on crack-tip-opening 
displacement (CTOD). Thus, Eq (17) is restricted to 

w - a - 0,53 ~ °‘ 65ct + °» 32ct2 (18) 

so that the influence of backface yielding on CTOD is less than 2 percent. The 
influence of backface yielding on plastic-zone size is not included in Eq (17). 

The CTOD (2V a ) for the compact specimen was calculated by adding the 
displacement at the tip of the physical crack length (a) due to the pin load 
and due to the uniform stress. Again, an equation was fitted to these results 
and was 


2V E 

a 


= 8 


2n F 1 F 2 " 


w a 

o 

irP 


✓8 


(19) 


where 


F 1 = (2 + A ) (o . 886 + 4.64A - 13.32A 2 + 14.72A 3 
- 5.6A 4 )/(1 - A) 3/2 

f 2 = 1 + b iY + b 2 y 2 


= -1.25 + 9.76A - 20.15A 2 + 16.62A 3 


B 2 = 0.64 - 4.34A + 10.24A 2 - 8.46A 3 

j[2 6 (l t fll + A 2 ) + ( 1 -M(5 + Ai-3A 2 )]| e 2 + (i - a)b 

+ (1 - A) 2 (3 - Aj + 3A 2 ) |jtn (/B + /B + 1 - A ) - £n /I - 


A) 


3/2 
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Aj = 3.57 + 12.5(1 - X) 8 

A 2 = 5.1 - 15.32X + 16.58X 2 - 5.97X 3 

H 2 = 1 + h^y + h 2 y 2 

hj = 0.666 + 0.796X 2 + 12.36X 4 

h 2 = 0.084 + 2.62A 2 - 14.08X 4 

6 = p/w, X = (a + p )/w, and y = p/(a + p). 

Equation (19) is within about ±1.5 percent of the current collocation 
results for 0.3 £ X £ 0.8, y £ 1 - 0.2/A, and p/(w - a) £ 0.5. Figure 8 
shows the normalized displacement against plastic-zone-size-to-ligament ratio 
for various a/w ratios. The curves show the results from Eq (19) and the 
symbols show results from Terada [9]. The results agreed well. The slight 
discrepancy at an a/w ratio of 0.3 was probably due to neglecting the pin- 
loaded holes in Terada 's analysis. Again the dashed curves indicate that back- 
face yielding is occurring, but its influence on CT0D is less than 2 percent if 
the plastic-zone size is less than that given by Eq (18). 

CONCLUDING REMARKS 

A two-dimensional elastic boundary-collocation analysis was used to ana- 
lyze the compact specimen subjected to various loading conditions. The effects 
of the pin holes were included. The loading cases were: pin-loaded holes, a 

pair of concentrated forces acting on the crack surfaces, and a uniform pres- 
sure acting over a segment of the crack surfaces. Stress-intensity factor (K) 
equations were developed for each of these loadings. Crack-surface displace- 
ment (V) equations were also developed for pin-loaded holes and for some 



restricted configurations of the uniform pressure case. The K-equations were 
within ±0.5 percent and the V-equations were within ±1 percent of the colloca- 
tion results. Some of these equations were used to develop algebraic equations 
for plastic-zone size and crack-tip-opening displacement (CTOD) for the Dugdale 
model. The plastic-zone size and CTOD equations were within 1 and 1.5 percent, 
respectively, of the collocation results. 

The results from this study on the compact specimen can be used in various 
analyses based on the Dugdale model concept. Some examples are the critical 
CTOD fracture analyses and the crack -closure models. 
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APPENDIX 

STRESS-INTENSITY FACTOR EQUATION FOR CONCENTRATED FORCES APPLIED 
TO CRACK SURFACES ON COMPACT SPECIMENS 

The stress-intensity factors for a compact specimen subjected to a pair of 
concentrated forces on the crack surfaces. Fig. 2(b), were obtained from a 
boundary-collocation analysis. An algebraic equation was fitted to the numeri- 
cal stress-intensity factors. 

For the compact specimen configuration, consider a semi -infinite crack 
along the x-axis in an infinite plate subjected to a pair of concentrated 
forces, Q, as shown in Fig. 9. The dashed lines Li and L 2 define the 
boundaries of the compact specimen. The Muskhelishvili [14] stress functions 
for this configuration are 

*(z) = $q(z) + ( z ) + *2^ 

*(z) = f Q (z) + ^(z) + t 2 (z) 

The subscripts denote functions which are needed to satisfy conditions for the 
concentrated forces on the crack surfaces and to approximately satisfy condi- 
tions on boundaries Lj and L 2 , respectively. 

The stress functions for a semi-infinite crack in an infinite plate sub- 
jected to a pair of concentrated forces on the crack surface were derived from 
equations given in Ref 15. They are 

v z > = v z > = 7 tan " 1 {F ( 2 D 

The stress functions used to approximately satisfy boundary conditions on 



the external boundary Li are 
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* x (z) 

^ x (z) 


= /z 


N 

I 

n=l 


n 


n 1 N 

2 1 A B » z 


n 


( 22 ) 


where the coefficients A n and B n are real. The stress functions will, 
of course, produce stresses on the internal boundary l_ 2 . 

The stress functions used to approximately satisfy boundary conditions on 
boundary L 2 are 



where C n , C n , D n , and D n are real. In these stress functions, 
poles 2 p| and zp, were located at the centers of the two holes (see 
Fig. 9). The stress functions in Eqs (22) and (23) automatically satisfy the 
conditions of stress-free crack surfaces and the single-valuedness of displace- 
ment condition for multiply connected regions. The conditions on boundaries 
Ll and L 2 were approximately satisfied by the series solution using the 
boundary-collocation method described in Ref 10. From a convergence study, N 
was selected as 40 and M was 20. 

The stress-intensity factor is given by 


K = 2/Zv lim /z $' (z) 
z+0 


(24) 
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They were calculated for a wide range in configuration parameters. An empiri- 
cal equation was fitted to the numerical stress-intensity factors and is 


K = 



(25) 


The equation for the correction factor, G, is 
G(o,a) = (l + A x A + V 2 )/(l - a) 3/2 

= 3.57 + 12.5(1 - a) 8 


(26) 


A 2 = 5.1 - 15.32a + 16.58a 2 - 5.97a 3 


where a = a/w, A = b/w', 0 £ A £ a, and 0.2 £ a £ 0.8. Equation (25) is 

within ±0.5 percent of the boundary-collocation results for 0.3 £ a < 0.8 and 
within ±1 percent for 0.2 £ a < 0.3. For clarity, only the results for a/w 

ratios of 0.3, 0.5, and 0.8 are shown in Fig. 10. The curves were calculated 
from Eq (26). 
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(a) "Dugdale" yield model. (b) Equivalent elastic loading. 

Fig, 1 Dugdale model for the compact specimen. 
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(a) Pin-loaded holes. (b) Concentrated forces. Cc) Uniform pressure. 


Fig. 2— Compact specimen configuration subjected to various loading. 
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Fig. 5— Normalized crack-surface displacement for compact specimen with 
concentrated forces on crack surface at b/a = 0.5 for various 
a/w ratios. 
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Fig. 6— Normalized crack-surface 
uniform pressure applied 
th/a ratios. 
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Fig. 7— Comparison of normalized plastic-zone size against 
applied load for various a/w ratios. 
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Fig, 8 Comparison of crack-tip-opening displacement against 
plastic-zone size for various a/w ratios. 








Fig, 10— Boundary-correction factors for compact specimen with pai 
of concentrated forces, Q, applied to crack surface. 
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